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Abstract 



O 
CT 

The goal of this paper is to propose an approach to the formulation of dynamics for 

causal sets and coupled matter fields. We start from the continuum version of the 

action for a Klein-Gordon field coupled to gravity, and rewrite it first using quantities 

that have a direct correspondent in the case of a causal set, namely volumes, causal 

relations, and timelike lengths, as variables to describe the geometry. In this step, 

the local Lagrangian density C(f;x) for a set of fields / is recast into a quasilocal 

expression £o(f;p,q) that depends on pairs of causally related points p ~< q and is a 

function of the values of / in the Alexandrov set defined by those points, and whose 

limit as p and q approach a common point x is C(f;x). We then describe how to 

OO . discretize Co(f;p,q), and use it to define a causal-set-based action. 

O 



1. Introduction 

A causal set pQ (for recent reviews, see Refs [21 OH]) is a locally finite partially ordered set, 
i.e., a set S with a partial order relation -< in which for every two elements p and q, the 
interval, or Alexandrov set, a(j>, q) := {r | p -< r -< q} is finite. The causal set approach to 
quantum gravity uses causal sets as the only fundamental structure for the description of the 
geometry of spacetime, and is based on the observation that, if we interpret the partial order 
as the set of causal relations and define the volume of a subset of a causal set to be simply 
the number of elements in it, we obtain a discrete version of all variables needed to describe 
a Lorentzian geometry. In the context of the continuum, it is relatively simple to see that a 
Lorentzian metric tensor field is equivalent to a light-cone structure and a volume element at 
every point, and in a causally well-behaved Lorentzian manifold the equivalence carries over 
to the causal structure. More precisely, the topology, differentiable structure, and conformal 
metric of any past and future-distinguishing spacetime are uniquely determined by the partial 
order defined by the causal relations among its points [3 a similar result holds when a 
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countable, dense set of randomly chosen points in the manifold is used [TJ. It is conjectured 
that, if the causal relations are known among points randomly chosen with finite density p, 
then the conformal structure can be approximately determined, up to variations on small 
scales compared to the point density, but in addition the volumes of spacetime regions (and 
thus the whole metric, up to a global scale factor) can be determined by counting points in 
them, with the same type of approximation. These kinematical considerations tell us that, 
although causal sets are not geometrical objects, and most of them do not correspond to 
Lorentzian geometries in any natural way, those that are "manifoldlike" in an appropriate 
sense are likely to define an effectively unique large-scale geometry. 

Formulating a reasonable dynamics for causal sets has remained an elusive goal. Various 
proposals have been made, of which the most studied one is the classical sequential growth 
dynamics of Rideout and Sorkin [8], in which causal sets ("causets" for short) are grown 
one element at a time, with each new element appearing to the future of some subset of the 
existing set, the specific probabilities for the choice of that subset being characteristic of the 
model used. Thus, each model gives a probability distribution on the set of iV-element causets 
after N elements have been generated. In this approach, if one imposes suitably formulated 
discrete versions of general covariance and causality, each possible model is characterized 
by a sequence of parameters, which allows the whole set of models to be systematically 
studied. As a result, some properties of causal sets that arise from these models are known, 
in particular for the special type called "transitive percolation" (or "random graph orders" 
in the mathematics literature) [9], but also for more general models [10] and, although any 
causal set can be generated by almost any one of the models, there are indications that the 
probability that a manifoldlike causal set emerges goes to zero as N — > oo [llj. 

In this paper, we will follow a different approach to the dynamics of causal sets, which 
is in principle more straightforward: We seek a way to write down the usual continuum 
action for gravity, and possibly matter fields, where the metric has been replaced by causal 
relations and volumes and the additional fields expressed in terms that have a well-defined 
causal set correspondence, and then discretize the resulting expression. In practice, this task 
is not as simple as it sounds. First of all, while using the volume element as a variable can be 
thought of as simply identifying the occurrence of y/—g in the action, it would not be useful 
to think of the locally defined conformal structure, \g\~ 1 ^ d g^u m d spacetime dimensions, 
which is tensorial rather than combinatorial in nature, as the variable representing the causal 
structure, for the purpose of discretizing the theory. Instead, we will use as variable the actual 
partial order p -< q, or q G J + {p), even in the continuum. This, however, will introduce some 
non-locality in the theory; we will elaborate on this point below. The difficulty involved in 
the inclusion of matter fields in this framework depends on the type of field considered, but 
the case of a scalar field is simple enough that we shall treat it first, even before the case 
of gravity itself; results for other fields will be discussed in separate articles [12]. The final 
issue related to our approach consists in specifying how to use the action we obtain in a 
dynamical principle. Although fully addressing this issue is beyond the scope of this paper, 
we will include some remarks on it. 
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2. Our Program 

To illustrate the issues that arise when writing a form of the action for a set of fields that 
is meaningful in the causal set context, consider the example of a Klein-Gordon field <b 
minimally coupled to a metric g^ v (considered as fixed, for now), with Lagrangian densitjo 

£ KG (^, <M = f Gr d,ct> d v 4> - m 2 <p 2 ) . (1) 

If our goal was simply to discretize the continuum theory defined by Eq (pQ), for example as a 
computational tool, we could introduce a regular lattice and convert the derivatives into finite 
differences using standard methods. In a fundamentally discrete theory this would not be a 
reasonable approach, because a regular lattice introduces extra structure in the theory and 
breaks its covariance. The only way to overcome this difficulty with the discretization process 
is to use a "random lattice" [13J; we will in fact view a causal set as a random discretization of 
a Lorentzian manifold (M, g^ u ) for some purposes. On a random lattice, however, covariant 
versions of finite-difference Lagrangian densities typically lead to divergent expressions. For 
example, we could try to replace the derivatives in Eq (JTJ) with a sum of terms of the type 
(0(9) ~~ 0(p)) 2 /( r (P; ?)) 2 ; where r(p,q) is the length of the longest future-pointing timelike 
curve from p to q, over all timelike "nearest neighbors" q of p (spacelike distances would 
not carry over well to the causal set context). But in a random lattice in a Lorentzian 
spacetime there are no "nearest neighbors", in general. Even in Minkowski space of any 
dimensionality d > 2, for any p G M. and e > the set of points q for which r(p, q) < e 
is a region of infinite volumep so there are infinitely many random lattice points within 
an arbitrarily small Lorentzian distance, and none of them minimizes that distance. Thus, 
since the terms we would be adding are all non-negative, any non-trivial, Lorentz-invariant 
definition of this type would give divergent results, and the divergence cannot be removed 
by some renormalization procedure without violating Lorentz invariance. 

The way to avoid the above infinity in £kg is to restrict ourselves to summing the terms 
(0(9) — 0G°)) 2 /( r (P> l)) 2 only over pairs of points within some finite region around x. In the 
context of causal sets, the simplest kind of finite region we can define using the partial order 
is an Alexandrov set. Thus, we will try to reconstruct g^d^cj) d u <p(x) by solely looking at the 
behavior of <ft i n the interior of an Alexandrov set centered at x. It may at first appear that 
this cannot work, because choosing one particular Alexandrov set is equivalent to choosing 
a preferred frame at x, the one in which the axis of the Alexandrov set coincides with the t 
direction, and our results will therefore not be covariant. However, what we are trying to do 
is determine the value of a scalar by using values of quantities defined in a particular frame, 
which is certainly possible. In fact, as we show in Sec 2, if is a different iable function, then 
it is possible to find a small enough Alexandrov set a(p, q) such that 

r dM) „ A (M^M L + _£_ iS m _ mf , (2) 

1 Although, strictly speaking, a Lagrangian density should include the volume element y/—g, in this paper 
we prefer to include the latter in the integration that gives the action S — J M d 4 x \f-~g C. 

2 It is the region between the future light cone of p, given by t — r if p is the origin, and the hyperboloid 
t = \lx 2 + e 2 . Since the surface area of a (d — 2)-dimcnsional sphere is r d ~ 2 for some positive constant 
Cd-2, the volume of the region is V — J^°(Vr 2 + e 2 — r) Cd-2 rd ~ 2 dr = (| g2 c d-2 r d ~ z + 0(e 4 )) dr = 00. 
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where A and B are appropriately chosen coefficients; a corresponding expression for gravity 
is derived in Sec 3. This result can be seen as analogous to finding an expression for a scalar 
invariant of the Faraday tensor in electromagnetic theory using a combination of the 
electric and magnetic fields in some reference frame, with appropriately chosen coefficients 
to make the result covariant. 

There is one important difference between the electromagnetic example and our situa- 
tion, however, which is that the quantities appearing in (J2J) do not have simple transformation 
laws under changes in a(p, q), except in the limit of infinitesimal sets, in which the A(f> change 
linearly under Lorentz transformations. As a consequence, the right-hand side is only covari- 
ant, and in fact only approaches the left-hand side, in the limit as p, q — > x, or a(p, q) — > {x}. 
Notice that this is not the same as the r(p,x), r(x,q) — > limit. The points p and q can 
be arbitrarily close to x, or to each other, in the latter sense and still have coordinates that 
differ by arbitrarily large amounts in any given reference frame, just by being sufficiently 
close to each other's light cone. Therefore, even smooth functions may not have well-defined 
limits in the sense of the Lorentzian distance. 

Once we have found an Alexandrov-set-based expression such as that in ([2]) for a La- 
grangian density that converges to the right value as p, q — > x, the remaining steps are: (i) 
Choose a specific pair (p, q) for each point x\ (ii) Write down a discretized version of the 
Lagrangian ([2]) where x is an element of a causal set that may or may not be manifoldlike; 
and (iii) Sum the contributions from all x to the action. The tricky thing about the first 
step is how to achieve the desired result in a covariant way; in Sec 4 we will describe one 
way of doing this, which consists in choosing the interval in which the field varies more 
slowly, in an appropriate sense. Part of the second step can be taken care of in an obvious 
way: we consider first the case of a causal that is manifoldlike, and embedded uniformly in 
a Lorentzian manifold with density p, in which case we can replace generic manifold points 
by embedded causal set elements, and integrals over Alexandrov sets by sums over elements 
contained in the corresponding intervals, each element with a weight p -1 ; since the resulting 
expression only contains quantities that are meaningful in a general causal set, we simply 
use the same one in the non-manifoldlike case. What is not so obvious in the causal set case 
is what to replace the p, q — > x limit by. That limit cannot be taken in the same way as in 
the continuum for the simple reason that the set is locally finite, and we restrict ourselves 
to using intervals with at least a certain minimum number N of elements, to minimize the 
effects of statistical fluctuations. The third step is then essentially trivial. 

3. Scalar Field 

Consider a massive Klein- Gordon scalar field in a curved (^-dimensional spacetime, d > 2, 
with Lagrangian density written in the usual form ([1]). We wish to derive an alternative 
expression for Ckg(9hv, (ft, x) in which the dependence on g^ v is expressed only through the 
causal structure and volume element. The action for gravity itself will be treated in the next 
section, but here we will lay some of the groundwork for that discussion as well. 

We are interested in calculating integrals of expressions involving (f>(x) over small Alexan- 
drov sets a(p,q), which in fact will approach a point in the continuum. We can therefore 
assume that a(p, q) is a convex normal neighborhood of its midpoint and expand all tensor 
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fields inside a(p,q) in Riemann normal coordinates x M based at that point, with the vector 
d/dx°\o pointing along the geodesic from p to q. For the scalar field we simply expand 

<j>{x) -<f>(p) = d fl <f)\ x» + o{x) , (3) 

for small x^, where "little o" is defined as usual by the property that o(x k ) vanishes faster 
than k-th powers of the coordinates as i* 1 — > 0. We can similarly expand all geometrical 
quantities. The metric itself, in Riemann normal coordinates, can be written in the form 

9„u — Vims + I ^W(°) xPx ° + °( x2 ) > ( 4 ) 

whose volume element is 

dV(x) = \g{x)\ 1/2 d d x = (l + 0(x 2 ))d d x , (5) 

and we can define an h so that the coordinates of the endpoints p and q are given by 
(Th, 0, 0) up to terms of order h 2 or higher; specifically, we choose 

h:=lr(p,q). (6) 
What we wish to do is calculate the integral on the right-hand side of (J2J), 

(0(x)-0(p)) 2 d^ = (9 M 0^0)L [ {x»-p»){x v -p v )d d x + o{h d+2 ), (7) 



where a is the Alexandrov set of the points (=Fr/2,0, ...,0) with respect to the flat metric 
r]^ u , in terms of derivatives of <fi. Since a is just a pair of cones of height h in Minkowski 
space over a (d — l)-dimensional ball of radius h and volume 

2ir d / 2 

V(B d ~ 1 , h) = c d _i h d ~ l , with c d 



dT(d/2) v ' 

(for example, c\ = 2, C2 = it, C3 = 1 7r, c 4 = ~ 7r 2 , c 5 = ^ 7r 2 , ...), the volume of a is simply 

V :=J a d d x = 2- C -^ = k d T d , k d :=^L. (9) 

Expanding the product in the integrand on the right-hand side of ((Tj), and noticing that 
integrals of odd functions over symmetric domains such as cto vanish, we can rewrite (J7J) as 

(<j>(x) - <j>(p)) 2 dV x = J2(d^\o) 2 ( [ (xn 2 d d x+(p») 2 [ d d x 



= fee + 1 M W) 2 + 7 d)1 r d + 2 (d fc 0«9 fc 0) , (10) 
where a summation over k = 1, d — 1 is implied, and we have defined the integrals 
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with 

J d := / (cos9) d d9 (13) 

J—k/2 

(thus, J\ = 2, J 2 = \ 7r, J3 = |, J4 = §7r, J 5 = i|, ...). To calculate the value of I^o, we 
can use the fact that the Alexandrov set «o is symmetric under 1 1— > — i, and each of the two 
cones it consists of can be sliced into t = constant hypersurfaces that are simply (d — l)-balls 
of radius h — \t\. To obtain the value of J^i, we further slice each ball into x\ = constant 
discs of radius \fjh — |^o|) 2 — (^i) 2 - 

We can now use Eqs (j2D , (EJ) and ( |T0i) to express the derivatives that enter £kg in terms 
of Alexandrov-set-based quantities. For the time derivative, (j3J) gives us 

do(j) = m - m (14) 

r 

(here and in the following we are neglecting higher-order terms in t _1 ), and for the spatial 
derivatives, substituting this last result into (TTUI) . we get that 
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dk^ = —^ (0(x)-0(p)) 2 d d x- 1 4 ^ ™» . (15) 

^,1 Ja(p,q) h,l ' 

These last two equations allow us to write the kinetic term in £kg as 

gT d^cj) d v <p = d (f) d (j) - d k <p d k <p 

J d,i J r 2 I dil T d + 2 J aM 

On the other hand, if we choose to use the volumes of the Alexandrov set instead of r to 
represent the separation between p and q, then gives us 



-) 



1/d 



(17) 



which, upon substitution into (JIB]) gives 

2/d 



™M-(i + *£*)(0 



(0(g) -m) 
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1 fv\ 1+2/d 

' 1 ' {(j>{x) -<j){p)) 2 d d x 

'a(p,q) 



IdA \h 



Thus, in four dimensions (for concreteness, in which case = vr/24, = 7r/ (2 -3-5), 
h,i = 7r /(2 4 ■ 3 2 • 5)), the Alexandrov set version of the scalar field Lagrangian density (JTJ is 

£kgH, \g\,<f>',P,q) (19) 
(0(g) 360 /■ (0(x) _ 0(p))2d 4 x+ i m202(p) 



r 2 (p,g) KT 6 (p,q) J a{p q) 
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up to o(l) terms, where, in the continuum context, we view the proper time as a function of 
the causal structure and volume element, obtained from Eq <^j, 

r( M ) = 2^<M>r. (20) 



Finally, the scalar field action in terms of 0, Alexandrov sets, and the volume element 
is obtained by taking the local limit and integrating over manifold points, 

£kgH, \g\A;x) = hm £kgK M,</>;p,<?) (21) 

Pi y * 

Skg= [ d d x\g\ x / 2 C KG {^\g\^;x) . (22) 
JM 

Notice that the quasilocal expression ffT9l) . viewed as a function of (-<, |<?|,t), without us- 
ing (I20T) to eliminate r, would in a sense resemble more the original C^Q^g^.fp.x), with 
derivatives expressed by terms containing lengths rather than volumes. However, since the 
continuum dynamics uses the local Lagrangian density, and in the limit fl2~Tl) the relationship 
f l2"U|) always holds, the set of variables (-<, \g\,r) is a redundant one, and this point of view 
makes no practical difference, even for the quantum theory. But in the discretized theory 
that is no longer true. Even for manifoldlike causal sets the relationship between r(p, q) and 
the length £(p, q) of the longest chain between p and q is not straightforward. Even for small 
distances, when those two quantities are approximately proportional in the mean, 

r(p,q)=K d (£(p,q)) , (23) 

up to curvature corrections, the coefficient K d is dimension- dependent and not known ana- 
lytically [IB] . Most importantly, the action needs to be extended to non-manifoldlike causal 
sets, for which £(p, q) as a variable is independent of the cardinality of the interval between 
p and q. 



4. Gravity 

We now wish to write down a similar form for the Einstein-Hilbert action for gravity, 

5ehW= / #cH4, (24) 

JM 



again in d > 2 spacetime dimensions. From Refs [Ml [15] we know that the generalization of 
the expression (Q for the volume of an Alexandrov set a of height r to curved spacetime is 

V(t) = k d r d [1 + (a d R + b d R 00 ) t 2 + o(r 2 )] , (25) 

where k d was defined in 0, the curvature is evaluated at the midpoint of a, 

d d 
ad = 24(d+l) (d + 2) ' d = 24(d+l) ' { } 
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and the o(r 2 ) terms contain higher powers of the curvature tensor. This equation can be 
rewritten as 

Yill-l = (a d R + b d Roo) r 2 + o(r 2 ) , (27) 

where the left-hand side contains the kinds of variables we want to use in the causal set 
context, while the right-hand side contains R, the Lagrangian density, and Rqq. Therefore, 
the idea is to use (|27[) to obtain two equations involving the same curvature components, 
and solve them for R. This will involve applying (|27|) to different Alexandrov sets expressed 
in the same coordinates, so we start by writing it in an equivalent covariant form, in terms 
of the endpoints p and q of a, 



v(p, g) 

k d r d (p,q) 



1 = {a d Rg, u + b d R pv ) (<f - p") (<f - p v ) + o(r 2 ) . 



(28) 



One of our two equations is (12 7p itself, and we obtain the other one by writing it for an 
Alexandrov set a(p,x) and integrating over all x G a(p,q). To leading order, this gives 



(p,q) \k d T d {p,x) 



1 J d d x = (a d Rr]^ v + b d R^) / ( 

/ J a(p,q) 



(x" - jf) (x u - p v ) d d x (29) 



a d R I [(x°) 2 -x k x k 

' a{p,q) 



+ p^] d d x + b d [ [Ro (x° - p°f + Ru {x l f] d 

{a d R [/ di0 - (d - 1) J d)1 + \ k d j + b d [Roo {hfl + \k d ) + Ru J d ,i] } r d+2 . 
^R a d (^I dfi - (d - 1) I dt i + \ k d ^ - b d I dt i + .Roo b d (^I dfl + h,i + \ k d ) } r 



-d+2 



where in the last step we have used J2i^-a = ~R + -^oo- The two Eqs (|27l) and (J29l) form 
a linear system for the unknowns R and i?oo- The determinant of the coefficients of those 
variables is 

D=(a d d + b d I d , 1 )r d+2 , (30) 
and when we solve the system for R = £eh we obtain the quasilocal expression 



£eh " D 



V{j) 
k d r d 



l) (l d ,o + Id,i + lk d )r d - ! f 



k d r d (p,x) 



- 1 )d d x > . 



(31) 



An important difference between this expression and the corresponding one for the scalar 
field, Ckg, is that now V(p, q) and r(p, q) are not equivalent variables, and therefore we 
cannot consider (-<, |g|, r) a redundant set as we did in the scalar field Lagrangian. Therefore, 
strictly speaking, the Lagrangian £ EH in (}3~Tj) is not expressed in terms of volumes and causal 
relations alone, since we need timelike lengths to evaluate it. While it is possible that an 
expression purely in terms of (-<, \g\) can be found, in this work we will use the above £eh, 
since in the causal set context there is a way to estimate timelike lengths independently of 
Alexandrov set volumes. 
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5. Discretization of the Lagrangian 

In the previous two sections we showed how to rewrite the action for a Klein-Gordon scalar 
field and for general relativity using just the causal structure -<, the volume element \g\, and 
the Lorentzian distance r as geometric variables. The general form of the action was that of 
an integral over Ai of a Lagrangian density for a set of fields / (which include the matter 
fields as well as the geometrical variables) obtained as a limit C(f;x) = \im p ^ x £ (f;p,q) 
of an expression that depends on / inside an Alexandrov set a(p, q) centered at x. The 
quasilocal Lagrangian density C (f;p,q) can be readily discretized, as we will show in this 
section; in the next section we will propose an approach to defining a discrete version of the 
limiting procedure needed to build the action from the quasilocal Lagrangian density. 

If we assume that the relationship between a continuum geometry {/A, g pv ) and a causal 
set C, when the latter is "manifoldlike", is that C appears as a uniformly random set of points 
embedded in (A4,g^ u ) with density p, then the discrete version C(f;p,q) can be obtained 
simply by taking p and q to be embedded causal set elements with p -< q, replacing the 
Alexandrov set a(p,q) in M. by the interval between p and q in C, and converting every 
integral over an Alexandrov set into a sum over interval elements, each carrying a weight 
p _1 , with any matter fields being evaluated again using the embedded points. The specific 
way in which this last step is carried out depends on the type of fields under consideration, 
but in the case of a scalar field it is obvious, and for a Klein-Gordon field Eq (|T6l) gives 

AcgH, T,(/>',p,q) (32) 

= (l + jjg + jM Mg) ~ Pl V U {x) _ m ? - I m 2 Sip? 

V hi ) 2r»(p,g) 2/ dil r^?) p t4 2^m- 

where d and p are considered parameters, and we obtain different versions depending on 
whether we express r(p, q) in terms of the cardinality V(p, q) of the interval a(p, q) using 
( 1201) or in terms of the length £(p, q) of the longest chain between p and q using Eq ( |23l ). 



We can similarly discretize the gravitational field Alexandrov-set-dependent Lagrangian 
density, and for the Einstein-Hilbert action Eq (I3TT) gives 



C ER (^,T;p,q) 

V(p,x) 



k d r d (p,x) 



(33) 



where d and p are still considered parameters, but in this case we need to think of r(p, q) 
as expressed in terms of £(p, q) using r = n^i, the relationship that holds in a manifoldlike 
causal set. 

From general facts about causal sets, we know that the manifoldlike ones — the ones 
"on the manifold shell", so to speak — are a tiny fraction of the total, but they are the 
only ones we know how to use for the purpose of checking the correctness of the dynamics, 
so we will just extend the expressions just obtained for the manifoldlike case to all causal 
sets; this is similar to what we do, successfully, in the path integral formulation of particle 
mechanics and flat spacetime field theory, where most of the histories that contribute are 
nowhere differentiable and nothing like the ones we consider classically. 
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6. The discrete action 

We will now address the issue of selecting the "right" Alexandrov set around each point, for 
the purpose of summing over points and computing our actions. There are two difficulties 
here: (1) If we select an Alexandrov set with too few points, then random fluctuations will 
no longer make our approximation reliable, and (2) If the axis of Alexandrov set is too 
close to the light cone, then despite the fact that the Lorentzian distance between the two 
endpoints of Alexandrov set might be small, the difference between the coordinate values of 
these points may still be arbitrarily large, hence arbitrarily large fluctuations of the fields 
are possible. The way to address the first issue is simple: we have to set the condition 
that the number of elements in the Alexandrov set be greater than some fixed number, or 
equivalently, V(a(p,q)) > V\. 

As far as the second condition is concerned, this one is a little bit more tricky, because 
the notion of being "sufficiently far away from the boundaries of the light cone" is not 
relativistically covariant. But then the question is this: why is the fact that the notion is 
not covariant not a problem in practice when talking about slowly-varying fields? To try 
to answer this, we can imagine a situation where things don't work near the light cone of a 
point p. For example, a planet that is a million kilometers away from p moves parallel to 
the t axis, and eventually crosses the light cone of p. Consider a point q, that lies on the 
path of that planet and is not on the light cone of p but to its future, so close to it that 
its Lorentzian distance to p is very small. We can now make a Lorentz boost so that p and 
q both lie on the new t axis. In this case, the picture we will see is that q is close to p 
coordinatewise, but this didn't help our situation because the planet was moving so fast in 
the new coordinates that despite the fact that very little time had elapsed between p and q, 
it still "had time" to hit q despite having been very far away from p. 

Thus, what we see is that the question we need to address is really not one about 
relativity and covariance, but rather how come the scenario just described is not a problem. 
In fact, the same issue arises in Newtonian physics and Galilean invariance, the only difference 
being that in that case the allowed range of velocities is infinite. In both cases, it would 
seem that we should run into the problem that, since the invariance group is not compact, 
if a planet were to pick a random velocity, it would pick with probability 1 one outside any 
given range of velocities. In special relativity, this means that if things were truly random, 
most planets would have been moving arbitrarily close to the speed of light, thus causing 
the disturbances we just described. So the fact that we see our spacetime as relatively 
continuous, without being disturbed by intruding planets, is an expression of the simple fact 
that in our universe, because of interactions, motion is not totally random, and the behavior 
of the matter in it selects preferred frames. What we have to do is find a way of selecting a 
preferred frame, in our case a preferred Alexandrov set, based on the behavior of fields near 
each point, and if we are able to do this covariantly we will not violate relativity. 

One way to select an Alexandrov set based on the behavior of the fields is to choose 
the one that has the smallest fluctuation in the Lagrangian density C, with the constraint 
already mentioned that its volume should be greater than some fixed V\. In order to do that, 
we have to be able to measure the fluctuations of C. If £ was a local function, then we could 
just find the point with the largest C and the one with the smallest £ in a given Alexandrov 



10 



set, and look at the difference between the two values. But our Lagrangian densities C are 
quasilocal functions of whole Alexandrov sets. Thus, instead of using points, we have to 
use smaller Alexandrov sets embedded into our large Alexandrov set. But then again we 
have to be careful: if the "small" Alexandrov sets are too small, then due to the fluctuation 
effects the results would no longer be reliable. So we have to again impose a lower bound 
on the volume of the "smaller" Alexandrov sets we are considering, which we will call V 2 . 
This means that the lower bound for the "large" Alexandrov set we are looking at would 
have to be "larger" than we had expected before: in addition to V\ being large enough for 
the "large" Alexandrov set to be reliable, it should also be large enough for there to exist 
V 2 <C Vi such that the Alexandrov sets of volume V 2 will also be large enough to be reliable. 

Let's make what we have just said more precise. Suppose we have a set of fields /, which 
includes both the gravitational field, -<, as well as any relevant matter fields. Our expression 
for the Lagrangian is a function both of / and an Alexandrov set a(p, q): 

C = C(f;a(p,q)). (34) 

The fluctuation of C over the Alexandrov set based at points p and q is then given by 

fluct(£, /; V 2 , a(p, q)) = max{£(/; a(r, s)) | r, s G a(p, q); V(a(r, s)) > V 2 } 

- mm{C(f;a(r,s)) | r, s G a(p,q); V(a(r,s)) > V 2 } . (35) 

For any point p, define the set of points q such that the Alexandrov set a(p, q) minimizes 
the above fluctuation, subject to the constraint that its volume be greater than V\. 

Q(C, f,p, Vl, V 2 ) = {qy p\Vr y p such that (V(a(p,r)) > V x , 

fluct(£, /; V 2 , a(p, r)) > fluct(£, /; V 2 , a(p, q))} . (36) 

Notice that in practice we expect Q to consist of a single point, in most cases. 

Finally, the Lagrangian density at a point p is given by the average of the values deter- 
mined by the set of all "smoothest" Alexandrov sets; in other words, we are averaging over 
all Alexandrov sets a(p, q) with q G Q(C, f,p, V±, V 2 ), 

rtt t, v \/\ - ^geQ(A/,P,vi,V2)^(./ , ' a (P^)) ,o 7 a 

Vl ' ' y " uWJ^VuW) • (37) 

where (!{•••} stands for the cardinality of a set. 



6. Conclusions 

In this paper we have obtained expressions for the action for a Klein-Gordon scalar field 
and for gravity in the causal set approach to quantum gravity. This was achieved in three 
steps. The first one, involving just the continuum theory for a set of fields /, consisted 
in finding a quasilocal Lagrangian density £ (f;p,q) defined for pairs of causally related 
points p -< q, such that its limit as p, q — > x is the desired local density C(f;x), and the 
metric appears in Co only through the causal structure -<, the volume element \g\, and the 
Lorentzian distance r. The second step was the discretization; this is logically the pivotal 



11 



point in the process, but in practice it turned out to be entirely straightforward, at least 
for the theories considered here. The key feature of our approach is that the discretization 
is carried out at the level of the quasilocal, Alexandrov set version £ of the Lagrangian 
density, because it is only there that a quantity governing the dynamics depends on the 
spacetime geometry only through the triple (-<, |<?|,t), each element of which has a causal 
set counterpart. The third step was the limiting procedure which gave us a discrete version 
of the local Lagrangian density, whose sum over causal set elements is the discrete action. 

The expressions we found do not constitute yet a full proposal for the discrete dynamics 
of causal sets and matter fields, for various reasons, in addition to the fact that we have not 
spelled out how they are to be used. For one, they depend on two integer numbers V\ and 
V2, for which we have not specified how to choose values. In addition to those parameters, 
they also include some ambiguities, including (i) the possibility of replacing r(p, q) (but not 
r(x, y)) in some or all of the places where it appears in (1521 and (1331) by V(p, q) using Eq 
(120]) . and (ii) different choices for how to extend the action off the "manifold shell". It is 
also possible that our expressions are simply an intermediate step in the derivation of other, 
more fundamental expressions in causal set terms. Finally, we should keep in mind that 
in the continuum one normally has to add boundary terms to the Einstein-Hilbert part of 
the action to satisfy certain criteria such as finiteness for asymptotically flat spacetimes and 
additivity, and those additional terms may change the structure of the causal set version of 
the action once they are taken into account. 

On the other hand, the fact that the dynamics of matter fields can be treated together 
with gravity along the same lines is a definite advantage of our approach. The approach raises 
interesting questions, such as whether the discrete action, obtained assuming the causal set 
was manifoldlike, will lead to such causal sets being dynamically preferred. The role played 
by Eq (1201) in our construction, as well as the fact that we expect the discrete Lagrangian 
density £ to approach a value independent of V\ and V2 as their values decrease (aside from 
"statistical fluctuations") also singles out one of the possible necessary conditions we can 
use to identify manifoldlike causal sets. If a causal set satisfies the property that, as two 
elements p -< q approach a common element x, the longest chain between them and the 
cardinality of the interval they define are related by Eq (1201) for some d that is independent 
of x (once statistical fluctuations are taken into account), we may not be able to conclude 
that the causal set is manifoldlike, but the different versions of the actions (|32|) and (1331) 
obtained from the continuum will coincide. 
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